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ABSTRACT 

We report the discovery of a substantial population of star-forming galaxies at 
3.0 ^ z ^ 3.5. These galaxies have been selected using color criteria sensitive to the 
presence of a Lyman continuum break superposed on an otherwise very blue far-UV 
continuum, and then confirmed with deep spectroscopy on the W. M. Keck telescope. 
The surface density of galaxies brighter than IZ = 2b with 3.0 ^ z ^ 3.5 is 0.4 it 0.07 
galaxies arcmin"^, approximately 1.3% of the deep counts at these magnitudes; this 
value applies both to "random" fields and to fields centered on known QSOs. The 
corresponding co-moving space density is approximately half that of luminous {L ^ L*) 
present-day galaxies. Our sample of z > 3 galaxies is large enough that we can begin 
to detail the spectroscopic characteristics of the population as a whole. The spectra of 
the z > 2> galaxies are remarkably similar to those of nearby star-forming galaxies, the 
dominant features being strong low-ionization interstellar absorption lines and high- 
ionization stellar lines, often with P-Cygni profiles characteristic of Wolf-Rayet and 
0-star winds. Lyman a emission is generally weak (< 20 A rest equivalent width) and 
is absent for >50% of the galaxies. We assign approximate mass scales to the galaxies 
using the strengths of the heavily saturated interstellar features and find that, if the line 
widths are dominated by gravitational motions within the galaxies, the implied velocity 
dispersions are 180 < cr < 320 km s-\ in the range expected for massive galaxies. The 
star formation rates, which can be measured directly from the far-UV continua, lie in the 
range 4 — 25 h'^^ Mq yr~^ (for go = 0.5), with 8.5/i^q^ Mq yr"-*^ being typical. Together 
with the morphological properties of the z > 3 galaxy population, which we discuss in a 
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companion paper (Giavalisco et al. 1996), all of these findings strongly suggest that we 
have identified the high-redshift counterparts of the spheroid component of present-day 
luminous galaxies. In any case, it is clear that massive galaxy formation was already 
well underway by z ~ 3.5. 

Subject headings: cosmology: observations-galaxies:formation-galaxies: evolution- galax- 
ies: distances and redshifts 

1. INTRODUCTION 

Several years ago we embarked on a program to assess the presence of "normal" star-forming 
galaxies at redshifts z > 3 using a custom suite of broad-band filters (dubbed Un, G, and TZ- 
see Stcidel & Hamilton 1993) designed to reach very faint apparent magnitudes with accurate 
colors. The technique rests on only two assumptions about the spectra of high-redshift galaxies: 
(i) a. roughly flat spectrum {f^, oc u^) in the far-UV where the flux is dominated by emission from 
massive stars and, (ii) & pronounced spectral "break" at 912 A (the Lyman limit in the galaxy rest 
frame) due to the combination of the intrinsic stellar energy distributions, the inherent opacity of 
the galaxies to Lyman continuum photons, and the effects of intervening absorbers (see Steidel & 
Hamilton 1992, 1993; Steidel, Pettini, & Hamilton 1995 [Papers I, II, HI] for an extensive discussion 
of the technique; see also Madau 1995). The filters are "tuned" to selecting objects with redshifts 
in the range 3.0 ^ z ^ 3.5, where the Un passband falls mostly shortward of the Lyman limit and 
the G band is not severely blanketed by the Lyman a forest — at higher redshifts this blanketing 
can make it very difficult to distinguish a spectral break from an intrinsically red spectral energy 
distribution. 

While earlier work (Guhathakurta et al. 1990) had shown that Lyman break galaxies do not 
dominate the galaxy counts at faint apparent magnitudes, these constraints were too weak to rule 
out the presence of a substantial population of such objects among the more prominent foreground 
blue galaxies. Our first searches were directed at the fields of QSOs with known Lyman limit 
absorption systems at z > 3; these could then be used as "templates" to guide one (eventually) to 
galaxies in the general field, or at least to test the feasibility of detecting normal galaxies at such 
extreme redshifts. In Paper III we discussed together the properties of both QSO absorbers and 
candidates in the general field, since QSO absorbers arc apparently drawn from the population of 
relatively luminous (but otherwise normal) field galaxies, at least at smaller redshifts (sec Steidel, 
Dickinson, & Persson 1994). We found the surface density to be the same within the errors in all 5 
fields studied, and we reported a surface density of "robust" candidates for Lyman break galaxies 
of 0.5 arcmin"^ to an apparent magnitude limit oiTZ = 25. 

Since the completion of Paper III, we have been working on extending the UnGTZ imaging 
technique to more random (i.e., non-QSO) fields and, most importantly, on following up the candi- 
dates with deep spectroscopy on the Keck telescope. In this paper we report the successful results 
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of our first attempts at spectroscopy of the Lyman break candidates, new UnGTZ imaging in (much 
larger) random fields which confirms the surface density estimates given in Paper III, and the first 
measurements of K-hwd magnitudes for the z > 2> population. 

2. SPECTROSCOPY 

The spectra were obtained during the nights of 1-2 October 1995 (UT) and 29-30 October 
1995 (UT) with the W.M. Keck Telescope and the Low Resolution Imaging Spectrograph (Oke et 
al. 1995). By using multi-object slit masks, each ~ 4' by 7' (and in one case long-slit spectroscopy), 
we recorded the spectra of 25 robust candidates for z > 3 galaxies in three fields: two are QSO fields, 
0000-263 and 0347-383, and the third is the Hawah deep survey field SSA 22 (Lilly et al. 1991; 
Cowie et al. 1994). These observations were made in far from optimum conditions, being either at 
large airmass (typically 1.8-2.0 for the southern fields) or in grey time (SSA 22); consequently, the 
spectra we obtained do not reflect the full capabilities of the LRIS+Keck combination. 

Candidates in the 0000—263 and 0347—383 fields were selected on the basis of the photometry 
in Papers I and III, augmented by deeper images of the former field acquired more recently. For 
the SSA 22 field we used photometry obtained with the Palomar 5m telescope in 1995 August and 
described in §4 below. As the number of available slitlets exceeded the number of robust candidates 
which could be observed with each mask, we also observed four "marginal" candidates — that is 
objects whose colors placed them at or near the boundaries of the color criteria (in either G — IZ 
or Un — G) outlined in Paper III — for the purpose of possibly improving our selection criteria in 
future. We summarize the results of the spectroscopy in Table 1, which includes relevant data for 
all of the candidates, both robust and marginal, which have been attempted to date. 

Of the 23 objects identified, 16 are galaxies with 3.01 < z < 3.43 and 3 are faint QSOs. The 
other 4 identified objects were considered to be marginal candidates, in 3 cases because of red G — TZ 
colors and significant Un band detections (all three turn out to be halo subdwarf stars). The 4th 
case, which showed evidence for a spectral break but which had a Un — G color slightly bluer than 
our robust candidates, is a galaxy at z = 2.780. Of the 8 robust candidates with undetermined 
redshifts, none is actually inconsistent with being at z > 3; in no case do we find a secure spectral 
feature at a lower redshift (the spectral range covered includes [OIIJA3727 up to z < 1.4). Many 
of these objects exhibit continuum discontinuities (~ 0.3 — 0.6 magnitudes) which are indicative of 
the onset of the Lyman a forest. 

To summarize, roughly 70% of the objects which we considered to be robust candidates have 
been confirmed to have 3.0 < z < 3.5, and the other ~30% are consistent with the same range of 
redshifts .0 



We note that an additional Lyman break candidate from our survey has recently been confirmed spectroscopically. 
In Paper III we proposed that an 7?. = 25.0 galaxy 2'.'9 from the line of sight to Q2233+1310 is the z = 3.15 damped 
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3. THE SPECTRA OF z > 3 GALAXIES 

In Figure 1 we have reproduced some examples of Keck spectra of z > 3 galaxies, chosen to 
illustrate the variety of features encountered. In each case we have included for comparison a recent 
HST spectrum of the central starburst region in the Wolf-Rayet galaxy NGC 4214 (Leitherer et 
al. 1996). 

Qualitatively, the similarity between the high-redshift galaxies and local examples of starbursts 
is striking. In each case the dominant characteristics of the far-UV spectrum are: (i) a flat contin- 
uum; (ii) weak or absent Lyman a emission; (in) prominent high- ionization stellar lines of He II, 
C IV, Si IV and N V and (iv) strong interstellar absorption lines due to low-ionization stages of 
C, O, Si and Al. These stellar and interstellar lines are the most distinctive spectral features in 
many of the z > 3 galaxies; this makes confirmation of objects at z > 2 with optical spectroscopy 
much more promising than might have been anticipated, even in the absence of strong ultraviolet 
emission lines. The continuum specific luminosity at 1500 A of a typical z > 3 galaxy in our sample 
is Li5oo = W^^h^Q erg s^^ A^-*^ (for qo = 0.5; 3 times greater for qo = 0.05); this is 500 — 1500 
times higher (depending on the value of qo) than the knot of star formation in NGC 4214 observed 
by Leitherer et al. (1996) and ~ 30— 100 times higher than the brightest such knots seen in nearby 
starburst galaxies. Thus, for a 'normal' IMF, the far-UV continuum we see in the z > 3 objects is 
produced by the equivalent of 2 — 6 x 10^ 07 stars. 

A detailed discussion of the spectra is beyond the scope of this paper; here we limit ourselves to 
a few preliminary considerations. All of the galaxies in the sample are consistent with unreddened 
models of young star forming galaxies (Bruzual Sz Chariot 1994). The expected dust-free G — TZ 
colors are in the range 0.3 — 0.9 after accounting for the blanketing of the G band by the Lyman 
alpha forest (see Madau 1995, Paper III), the same as the observed range in our sample. Modulo the 
uncertainties in models of young galaxies (see, e.g.. Chariot 1996), we estimate that the maximum 
reddening allowed in the G — TZ color is ~ 0.4 magnitudes. Taking this as an upper limit, and using 
the "extragalactic" reddening curve of Calzetti et al. (1994), the extinction at observed TZ (rest 
~ 1600 A) would be < 1.7 magnitudes, or a factor of < 5. The corresponding optical reddening in 
the rest frame of the galaxies would be E{B — V)< 0.3 magnitudes. 

Despite the apparent lack of a large amount of dust, Lyman a emission is always much weaker 
(usually by factors of more than 10) than the ionization-bound, dust-free expectations, given the 
production of ionizing photons by massive stars which we measure directly from the UV continuum. 
Reasons for the preferential extinction of Lyman a emission have been discussed by, e.g.. Chariot & 
Fall (1993) and Chen &; Neufeld (1994); our observations are entirely consistent with the same low 
(but non-zero) dust content inferred to be present in the high redshift damped Lyman a aborption 



Lyman a absorber, as this is the best candidate Lyman break object in the small field wc observed around the QSO. 
Djorgovski et al. (1996) have serendipitously obtained a spectrum of this candidate showing Lyman a emission at 
the predicted wavelength. 
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systems (e.g., Pei et al. 1991, Pettini et al. 1994). In the spectra which do show Lyman a emission 

the typical rest-frame equivalent width is Vl^Lyo = 3 — 20 A, but apparently most z > 3 galaxies 
have Lyman a emission lines weaker than these values, in striking similarity to nearby star-forming 
galaxies (e.g., Giavalisco, Koratkar, &; Calzetti 1996). For example, C23 0000-263 in Figure 1 has 
among the largest implied star formation rates in our sample (as discussed in §6 below), and yet 
has no measurable Lyman a emission. 

The Lyman a emission lines and number density that we measure for our galaxies are consistent 
with the essentially null results of all previous searches for high-z galaxies based on this spectral 
feature (e.g., Thompson, Djorgovski, h Trauger 1995; Lowenthal et al. 1995; Pettini et al. 1995). 
Given the star formation rates we derive, the space density of the Lyman break galaxies (§6) is 
also consistent with null results from near-IR surveys (e.g., Pahre & Djorgovski 1995; Mannucci 
et al. 1994) for redshifted Ha, [OIII], and [Oil] lines. While the near-IR surveys can reach star 
formation rates comparable to those we observe, the volumes surveyed to date are too small for a 
significant chance of detection. 

While there is a great deal of variety in the strengths of the lines which are predominantly 
of stellar origin (C IV A1549, Si IV AA1393,1402 and He II A1640), we find these features to be 
generally weaker than in the spectra of present-day starbursts. This is probably an abundance 

effect; these lines are formed predominantly in the winds of massive stars, where both mass-loss 
rates and wind terminal velocities are known to depend sensitively on metallicity (e.g. Walborn et 
al. 1995). 

The strongest interstellar lines indicated in Figure 1 have typical rest-frame equivalent widths 
IVq — 2 — 3.5 a. While the interstellar medium of these galaxies has obviously undergone some 
chemical enrichment, it is not possible to deduce metallicities from our spectra, since the absorp- 
tion lines in question are undoubtedly heavily saturated. Under these circumstances, values of 
metallicity anywhere between 1/1000 of solar and solar are compatible with the line strengths and 
higher-resolution observations of intrinsically weaker lines are required to measure element abun- 
dances (Pettini &: Lipman 1995). The equivalent widths of saturated absorption lines are much 
more sensitive to the velocity dispersion of the gas than to its column density. The observed 
VFo ~ 2 — 3.5 A correspond to FWHM > 400 — 700 km s"-*^ which in turn imply approximate 
velocity dispersions a = 180 — 320 km s~^ (slightly higher values apply if rotation is the dominant 
effect). In any case, if these velocity spreads reflect primarily gravitationally induced motions in 
the large-scale interstellar medium of the galaxies, the masses implied are comparable to those of 
present-day luminous galaxies. Smaller masses would result if interstellar shocks, local to the star- 
forming regions, contribute to the line widths. Spectra of higher S/N and resolution are required 
to resolve this question. 
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4. THE SURFACE DENSITY OF LYMAN BREAK OBJECTS 

In parallel with this spectroscopic follow-up, we have been obtaining deep UnGTZ images of 
additional fields using the COSMIC camera at the prime focus of the Palomar 5 m telescope with 
two principal aims: (i) Improve the statistics on the surface density of Lyman break objects, by 
taking advantage of the much larger field of view available to COSMIC (9' 7 x 9' 7, compared to 
approximately 1'5 x l!5 of most of the observations reported in Paper III); and (ii) Extend this 
work to random fields — that is fields which do not include known QSOs — so as to assess if the 
surface density deduced in Paper III is typical of the general field. 

To this end we have chosen one of the Hawaii deep survey fields, SSA 22, which has also 
been observed extensively with post-refurbishment HST (e.g., Cowie, Hu, &: Songaila 1995a) in 
several pointings, all falling within our Palomar 5 m images. These new UnGTZ images reach depths 
comparable to those of the data in Paper III, and allow us to be complete in our identification of 
z > 3 galaxies down to 7^ = 25.0. We find a total of 31 robust candidates satisfying the selection 
criteria outlined in Paper III out of a total of 2340 objects detected to the same apparent magnitude 
level (the field size in the stacked images was 8' 8 x 9^0). We therefore deduce a surface density of 
Lyman break candidates of 0.40 ± 0.07 arcmin^^ (or 1.44 x 10^ per square degree), consistent with 
(and significantly more accurate than) the value ~ 0.5 arcmin"^ reported in Paper III. 

We conclude that: (i) Lyman break objects in the redshift range 3.0 < z < 3.5 represent about 
1.3% of all objects to 7^ = 25.0, and 2.0% of all objects in the magnitude range 23.5 < TZ < 25.0; 
and (ii) the density of z > 3 galaxies is not significantly higher in fields including bright high-z 
QSOs. 

5. THE NEAR IR EMISSION FROM z > 3 GALAXIES 

In 1995 October we obtained deep Kg band images of a small subset of the z > 3 candidates 
using the Near Infrared Camera (Mathews & Soifer 1994) on the W. M. Keck telescope. A total of 
5 candidates was observed, with typical integration times of 6000s; 4 of the candidates are among 
those confirmed spectroscopically (see Table 1). The measured K magnitudes (typically measuring 
rest-frame B or V of the galaxies) of the non-AGN candidates (0000—263 C09 appears to be an 
AGN) range from iC = 21.3-22.1, with 2.3<7e-Er<3.2[3. These colors are consistent with 
models of continuous star formation which could have begun as early as ^ 1 Gyr prior to the epoch 
at which we observe them and are redder than would be expected if we were seeing instantaneous 
"bursts" of star formation (Bruzual & Chariot 1993). A consequence of adopting the maximum 
amount of reddening allowed by the UV colors of the galaxies (§4) is a reddening in 7?. — of ~ 0.8 
magnitudes; this would lower the ages significantly and formally allow single burst models younger 



The optical passbands are on an "AB" -normalized system, whereas the K magnitudes are the standard system. 
A flat spectrum {fv =const.) source has TZ ~ K — 1.85 
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than a few times 10 years. We regard such short lifetimes as unlikely as they would imply that 
we are seeing large numbers of galaxies all bursting simultaneously. 

6. IMPLICATIONS 

The confirmation of a substantial population of luminous, star-forming galaxies at z > 3 has 
considerable implications for our understanding of galaxy formation and evolution. One of the 
advantages of searching for high— z galaxies in the optical is that one observes directly the far- 
UV continuum produced by early-type stars (our TZ bandpass samples the rest frame continuum 
at Aq ~ 1600 A). Consequently, in the absence of dust (see §5), relatively minor assumptions 
are necessary to deduce the formation rate of massive stars and the accompanying production of 
ionizing photons; this is not the case for measurements of [Oil] line luminosities which are subject 
to uncertainties as large as a factor of ~ 5 (see Gallagher, Bushouse, & Hunter 1989). For the 
purpose of estimating star formation rates from the observed UV luminosities, we have made use 
of the calculations by Leitherer, Robert, & Heckman (1995) which are based on ultraviolet libraries 
of massive star spectra coupled with an evolutionary synthesis code. We have used the "continuous 
star formation", rather than "single burst", models; we consider the former case more likely, for 
reasons given in §5. For a Salpeter initial mass function with an upper mass cut-off of 80 Mq,]^ a 
SFR = IMq yr~^ produces a luminosity at 1500 A (in the rest-frame) L1500 = 10^*^'^ erg s^^ A^^. 
At z = 3.25 (in the middle of the redshift range to which we are sensitive to the detection of Lyman 
break galaxies), an apparent TZ = 24.5 (on the AB system) corresponds to L1500 = 10^^'^ erg s^^ A^^ 
{Ho = 50 km s^^ Mpc^^) for qo = 0.5, and L1500 = lO'*^'^ erg s~^ A^^ for qo = 0.05 . The population 
of Lyman break galaxies we detect is in the range 25.5 ^TZ ^ 23.5; therefore, the implied star 
formation rates range from (4 — 25)h^Q Mq yr~^ {qo = 0.5) to (12 — 75)/i^g^ Mq yr~^ (go = 0.05), 
with the weighted average being S.bh'^o i'^^h^o) ^0 y^~^ for qq = 0.5 (0.05). 

Assuming that we have uniformly probed the redshift range 3.0 < z < 3.5 in our surveys (an 
assumption that is supported by the results of the spectroscopy), the co-moving density of the star- 
forming galaxies is at least 3.6 x lO'^h^Q (6.7 x W~^hlo) Mpc~^ for qo = 0.5 (0.05), or about 1/2 
(1/10) of the space density of present-day galaxies with L > L* for qo = 0.5(0.05) (Elhs et al. 1996). 
The total star formation rate per co-moving volume produced by the observed population at 2; > 3 
is then 3.1 x lO^'^/iso (1.8 x lO^^/iso) Mq yr^^ Mpc^^ for qo = 0.5 (0.05). These numbers must be 
viewed as strict lower limits on the total star formation rate at z > 3 — the fraction arising in only the 
most actively star-forming objects; taken at face value, the star formation density that we observe 
is only ~ 25% of the total star formation rate seen at the present epoch (Gallego et al. 1995). For 
comparison, the star formation rates (per object) and star formation densities recently reported by 



The parameters of the IMF are relevant for extrapolating the observed rate of formation of massive and short- 
lived stars to the total SFR; for our purposes the slope of the IMF is the important parameter, while the upper mass 
cut-off has only a minor effect — see Table 2 of Leitherer et al. 1995. 
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Cowie, Hu, &; Songaila (1995b) at z ^1 are approximately the same as we have observed at z ^3; 
however, as we detail in Giavalisco, Steidel, & Macchetto (1996), the morphology of the objects 
hosting the large star formation rates seems to be entirely different in the two redshift regimes — 
the great majority of the z > 3 objects have very compact morphologies, with half-light radii in 
HST images of 0.2 — 0.4" (1.5 — 3h'^Q kpc for go = 0.5), and in most cases are relatively "regular" 
in appearance. As pointed out in Paper III, it will be very interesting to observe how rapidly the 
surface density of Lyman break objects increases at fainter apparent magnitudes. 

In summary, we have discovered a population of star-forming galaxies with redshifts 3.0 < 
z < 3.5 using a color-selection technique whose efficiency is very high, allowing the first systematic 
study of the nature of galaxies at such large redshifts. The space density, star formation rates, 
morphologies and physical sizes (see Giavalisco et al. 1996), masses, and early epoch of the galaxy 
population that we have isolated all support the possiblity that we are seeing directly, for the 
first time, the formation of the spheroidal component in the progenitors of present-day luminous 
galaxies. The fact that we detect such a substantial population using a flat-UV spectrum selection 
criterion suggests that dust obscuration may not be an important limiting factor in searches for 
high-redshift galaxies. In any case, our results demonstrate beyond any doubt that massive galaxy 
formation was well underway hy z = 3.5. 

We would like to thank the staff of the W.M. Keck Observatory, without whom these obser- 
vations would not have been possible. CCS acknowledges support from the Sloan Foundation and 
from the NSF through grant AST-9457446. We thank W. Sargent, M. Fall, and A. Phillips for 
comments on an earlier draft of the manuscript. We are especially grateful to Claus Leitherer for 
providing us with the spectra of nearby star-forming galaxies prior to publication. 
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Table 1. Summary of Candidates Observed Spectroscopically 
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25.00 


0.52 


1.52 




2.780 


<0 


0000-263 C28 


24.90 


0.70 


>1.90 








SSA22 COl 


25.50 


0.25 


>1.93 








SSA22 C03 


25.37 


0.12 


>1.95 








SSA22 CIO 


24.86 


0.21 


>2.26 




3.375 


<0 


SSA22 Cll 


24.82 


0.41 


>2.18 








SSA22 C12 


24.78 


0.43 


>2.13 




3.201 


3 


SSA22 C14 


24.64 


0.26 


>2.31 




3.401 


<0 


SSA22 C16 


24.62 


0.46 


>2.10 








SSA22 C19 


24.58 


0.39 


>2.50 




3.019 


22 


SSA22 C20 


24.52 


0.27 


>2.55 




3.019 


<0 


SSA22 C24 


23.45 


0.38 


>3.01 








SSA22 DIO 


21.75 


0.52 


2.56 




3.083 


(QSO) 


SSA22 Dll"^ 


21.53 


1.13 


3.26 




star 




SSA22 D12 


20.98 


0.87 


3.82 




3.352 


(QSO) 


2233+131 Nl 


25.02 


0.15 


>2.38 


[3.0] 


3.151^ 





"Spectroscopic redshift, when secure. 

'^Rcst equivalent widtli of Lyman a emission line for galaxies, in A; values <0 indi- 
cate net absorption at the position of Lyman a 

Object "G2" (Giavalisco et al. 1994)- Keclj spectrum shows presence of C IV, 
He II, O VI, and N V emission, hence G2 is probably an AGN 

''Spectroscopic confirmation from Djorgovski et al. (1996) 
'^Marginal candidate on the basis otUn — G and G — 71 colors 
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Fig. la,b. Examples of Keck spectra of z > 3 galaxies. Below each observed spectrum, we have 
plotted the spectrum of a star forming "knot" in NGC 4214, a nearby starburst galaxy (Leitherer 
et al. 1996), after shifting the spectrum to the measured redshift of the high-z galaxy. Some of the 
spectral features seen in local star-forming regions (both stellar and interstellar) have been labeled, 
for comparison with the high redshift objects (see text for discussion). 
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Fig. lc,d . — (continued) 



